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Abstract 

Wormhole spacetimes may be responsible for the possible loss of quantum 
coherence and the introduction of additional fundamental quantum indeter¬ 
minancy of the values of constants of nature. As a system which is known to 
admit such classical wormhole solutions, Einstein-Yang-Mills (EYM) theory 
is revisited. Since the classical wormhole instanton solution in this theory has 
been studied extensively thus far, in the present work, “quantum wormholes” 
are explored. Namely in the context of a minisuperspace quantum cosmology 
model based on this EYM theory, “quantum wormhole”, defined as a state 
represented by a solution to the Wheeler-DeWitt equation satisfying an ap¬ 
propriate wormhole boundary condition, is discussed. Finally, it is proposed 
that the minisuperspace model based on this theory in the presence of the 
cosmological constant may serve as a simple yet interesting system displaying 
an overall picture of entire universe’s history from the deep quantum domain 

all the way to the classical domain. 
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I. Introduction 


If one identifies the spacetime metric as the relevant gravitational held subject to the 
quantization, the topology of spacetime is expected to huctuate on scales of the order of the 
Planck length Ip = M~^. And of all types of conceivable spacetime huctuations, major con¬ 
cern has been the “wormhole conhguration” which is an object that can be loosely dehned as 
the instanton which is a saddle point of the Euclidean action making dominant contribution 
to the topology changing transition amplitude. Then one of the most crucial effects the 
wormhole (or more generally, these spacetime huctuations) may have on low-energy physics 
could be the possible ehective loss of quantum coherence [1-3]. For example, one may specu¬ 
late the situation where “baby universes” are pinched oh and carry away information. Then 
this kind of stereotypical information loss can lead to an ehective loss of quantum coherence 
as viewed by the macroscopic observer who cannot measure the quantum state of the baby 
universes. This issue of the possible loss of quantum coherence in (mostly the semiclassical 
approximations of) quantum gravity stood at the center of debates some time ago [1-3]. 
In a more careful and concrete analysis, however, these arguments concerning the possible 
quantum incoherence should be taken with some caution. For example, there is an argument 
by Coleman [3] that in the context of “many universe” interpretation (i.e., the third quan¬ 
tization formalism), the quantum incoherence will not be observed, namely the quantum 
coherence will be restored. In fact, of all the possible effects of the fluctuations in spacetime 
topology on the low energy physics, the most provocative one that immediately attracted 
enormous excitement was the advocation initiated by Baum [1] and by Hawking [1] and then 
rehned later by Coleman [3] that the wormholes have ultimately an effect of turning all the 
constants of nature into “dynamical random variables”. Thus this Baum-Hawking-Coleman 
(BHC) mechanism leads to a striking conclusion that the effects of (particularly) wormholes 
introduces into the low energy physics a fundamental quantum indeterminancy of the values 
of the constants of nature which can be thought of as an additional degree of uncertainty 
over the usual uncertainty in quantum mechanics. Unlike the issue of the loss of quantum 
coherence discussed earlier, however, this BHC-mechanism does not simply imply the elimi- 
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nation of the classical predictability of nature by the effect of quantum gravity. For instance, 
the BHC-mechanism actually leads to the prediction that the most probable value of the 
fully-renormalized cosmological constant is zero, indeed in exact agreement with the obser¬ 
vation. For the detailed arguments involved in the BHC-mechanism particularly concerning 
the most probable value of the fully-renormalized cosmological constant, we refer the reader 
to the literature [1,3]. But it seems fair to mention that the formulation of wormhole physics 
particularly the one put forward by Coleman is not without some inherent flaws. First, log¬ 
ically the Coleman’s wormhole physics formulation had faced some severe criticisms such as 
“sliding of Newton’s constant problem” [7] and “large wormhole catastrophe” [8] which seem 
to be unfortunately quite generic. However, it seems that the most fundamental and crucial 
difficulty associated with Coleman’s formulation is the use of saddle point approximation 
to the Euclidean path integral for quantum gravity since, as is well-known, the Euclidean 
Einstein-Hilbert action is not bounded below [9]. 

Now, then, note that both issues discussed thus far, namely the loss of quantum coherence 
and the determination of probability distribution for constants of nature (which are now 
random variables), are clearly based on the assumption that there are wormhole instantons 
as saddle points of the Euclidean action of the theory under consideration. Therefore unless 
one can demonstrate that there are large class of theories comprised of gravity with or with¬ 
out matter which admit Planck-sized wormhole instantons as solutions to the classical held 
equations, the discussion above on interesting effects of wormholes on low energy physics 
will lose much of its meaning. Unfortunately, thus far only a handful of restricted classes 
of theories are known to possess classical wormhole instanton solutions and they include ; 
Einstein-Kalb-Ramond (KR) antisymmetric tensor theory [4], Einstein-Yang- Mills theory 
[5] and Einstein-complex scalar held theory in the presence of spontaneous symmetry break¬ 
ing [6]. The classical wormhole instanton solutions and the semiclassical analysis of their 
ehects on low energy physics in these theories had been thoroughly studied in the literature. 
Here in the present work, we revisit the Einstein-Yang-Mills (EYM) theory of Hosoya and 
Ogura [5] which is known to admit classical. Euclidean wormhole instanton solutions. One 
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may wonder why anybody should repeatedly go through a well-studied theory like this one. 
Although the “classical” wormhole instanton as a solution to the classical held equations 
and its effects on low energy physics had been studied extensively, almost no attempt has 
been made concerning the serious study of “quantum” wormholes in the same theory. Thus 
in the present work we attempt to explore the existence of quantum wormholes in EYM 
theory. Now we describe briefly the approach we shall employ to explore the physics of 
quantum wormholes in our theory. And to do so, it seems necessary to distinguish between 
the definition of “classical” wormhole and that of “quantum” wormhole. In the classical 
sense, wormholes are Euclidean metrics which are solutions to the Euclidean classical field 
equations representing spacetimes consisting of two asymptotically Euclidean regions joined 
by a narrow tube or throat. In the quantum regime, on the other hand, and particularly in 
the context of the canonical quantum cosmology, quantum wormholes may be identified with 
a state or an excitation represented by a solution to the Wheeler-DeWitt equation satisfying 
a certain boundary condition describing the wormhole configuration. An widely-accepted 
such “wormhole boundary condition” is the one advocated by Hawking and Page [10]. And 
it states that wormhole wave functions are supposed to behave in such a way that they 
are damped, say, exponentially for large 3-geometries {-s/h oo) and are regular in some 
suitable way when the 3-geometry collapses to zero {y/h 0). Thus we shall construct a 
minisuperspace quantum cosmology model possessing SO(4)-symmetry based on the EYM 
theory and attempt to solve associated Wheeler-DeWitt equation. It will be demonstrated 
that there indeed is a solution to the Wheeler-DeWitt equation whose asymptotic behaviors 
satisfy the “wormhole boundary condition” stated above and hence can be identified with a 
wormhole wave function, namely a universe wave function for quantum wormholes. 

This paper is organized as follows : In sect.II, we recapitulate classical wormhole instanton 
solutions and the semiclassical analysis of their effects on low energy physics in EYM theory. 
Sect.Ill will be devoted to the study of quantum wormholes in this theory employing the 
approach described above. Finally in sect.IV, we summarize the results of our study and 
discuss their physical implications. 
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II. Classical wormhole physics in EYM theory 

We begin by reviewing the classical wormhole solution in EYM theory and its implications 
on low energy physics in a rather different context from the one in the literature [5]. 
Consider a system consisting of Yang-Mills (YM) theory of SU(2) gauge field coupled to 
Einstein gravity in the presence of the cosmological constant. This Einstein-Yang-Mills 
theory is represented by the Euclidean action [5] 

Ieym = / d^x^[K - -^R + - / d^xVh^iK - A'„) (1) 

Jm Idtt ^ JdM oTT 

where is the field strength of the YM gauge field 21“ with a = 1,2,3 being the SU(2) 
group index and gc being the gauge coupling constant. We also added Gibbons-Hawking 
gravitational boundary term on dM with h being the metric induced on dM and K being 
the trace of the second fundamental form of DM. 

Then by extremizing the action above with respect to the metric g^y and the YM gauge 
field A“, one gets the following classical field equations respectively 

Rjxu — -9fxuR + ( 2 ) 

with 

JC ^ 

= 0 ( 3 ) 

where and henceforth we use the notations = (9^5““ -|- 

2l^ = A^{—iT'^) and F^,, = F^^{—iT^) with being the usual gauge-covariant 
derivative and T“’s are SU(2) group generators obeying the SU(2) Lie algebra [T“,T^] = 
ie^^bcj'c normalization Tr(T“T^) = (5“^/2. 

Now we look for a Euclidean, SO(4)-symmetric wormhole solution which is an instanton 
describing the nucleation of a, presumably, Planck-sized baby (spatially-closed ] k = +1) 
FRW-universe. To this end, we begin by taking SO(4)-symmetric ansatz for the Euclidean 
metric and YM gauge field. First, the spatially-closed FRW-metric, which has the SO(4)- 
symmetry, is given by 
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( 4 ) 


ds"^ = g^ydx^dx'' = 

= a^[iV^(r)(ir^ + a^(r)(T“®(j“] 

where iV(r) and a(r) are lapse function and scale factor respectively and = {2/3nMp) has 
been introduced for convenience. Throughout this work, we shall work with non-coordinate 
basis with indices A, B = 0, a {a = 1,2,3) rather than with coordinate basis with indices 
p, z/ = r, 9, (j), tjj where (6*, 0, ip) are Euler angles parametrizing the spatial section of the 
manifold which has the geometry of S^. Then the non-coordinate basis 1-forms can be read 
off as 

= {e° = (TA^(r)(ir, e“ = (Ta(r)(T“} (5) 

where {<7°'} {a = 1, 2, 3) form a basis on the 3-sphere as just mentioned satisfying the 
SU(2) “Maurer-Cartan” structure equation 

da^ = (6) 

and can be represented in terms of 3-Euler angles 0<6<7i,0<(j)<27i and 0 < '^ < dvr, 
parametrizing 

= cosipdO -f sinipsinOdp, 

= sin-ipdO — cosipsinOdp, (7) 

= dip + cos6d<p. 

For later use, we also write down the associated vierbein and its inverse using the dehnition, 
= ef^dx^ , = 5^ and = <5^ where x^ = (r, 6*, 0, ip). 

0 0 0 

^cosip ^sinip 0 

simp —costp g 

asinO asinO 

— siri'ipcosO cosipcosO 1 

asinO asinO a 

Next we turn to the choice of SO(4)-symmetric ansatz for the YM gauge held and its 
held strength. Note that the SU(2) group manifold is also just like it is the case for 
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the geometry of the spatial section of the spacetime manifold. Thus one may choose the 
left-invariant 1-form {a“} as the “common” basis for both group manifold and the spatial 
section of the spacetime manifold. And this indicates that there is now “mixing” between 
the group index in the YM held and the non-coordinate basis frame index since we choose 
to work with non-coordinate basis. In the present work, we restrict our interest only to the 
magnetically-charged solution, i.e., F^- = 0 {i = 6, 0, 0). The motivation for this restriction 
will become clearer later on in the discussion of quantum wormhole physics. Namely from 
= 0, F“} {i,j = 9, 0, 0), we can write down the SO(4)-symmetric ansatz for 
YM held strength in non-coordinate basis indices 


= eJ^elF 


-0 


fJ,l' 




K = 


0 . 


( 9 ) 


Note that the choice of the ansatz for K above is motivated by the choice of ansatz for the 
Kalb-Ramond antisymmetric tensor held strength, F[^^x = h{T)e^^\ in the work by Giddings 
and Strominger [4] who looked for SO(4)-symmetric wormhole instanton solutions as well. As 
a matter of fact, this is an essential common structure shared by the Einstein-antisymmetric 
tensor theory considered by Giddings and Strominger [4] and the present EYM theory hrst 
studied by Hosoya and Ogura [5] that allows the existence of classical wormhole solution. 
Also note that this kind of starting point of ours is different from that in the original work 
of Hosoya and Ogura [5] in which they begin by taking the SO(4)-symmetric ansatz for the 
YM gauge potential rather than its held strength as we are doing here. The motivation for 
taking the ansatz for the held strength will become clearer shortly. 

Now, we would like to determine the condition that the function /(r) appearing in F^f 
should meet and the form of the associated YM gauge potential Af by demanding only 
that they be consistent with the dehnition = dA" -f- A and obey the Bianchi 

identity DF'^ = dF°' — A A*^ = 0. And to do so, we choose to take the temporal gauge 

A“ = 0 (or equivalently Aq = 0 in non-coordinate basis). This gauge choice is indeed natural 
since the spacetime is homogeneous and isotropic thus depends only on r-coordinate, there 
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is no gauge freedom associated with the space-dependent gauge transformation. Now we 
begin with the dehnition of the YM held strength = dA^- + A A^. On the one 

hand, using the YM gauge potential expressed in non-coordinate basis (in = 0 gauge) 
= A^^^dx^ = A^e^ [b = 1, 2, 3), one can evaluate the above dehning equation for the held 
strength with the help of Maurer-Cartan structure equation da°‘ = A given earlier. 

On the other hand, from our choice of SO(4)-symmetric ansatz for the YM held strength, 
one has = \F^^dx^ A dx'^ = A e'^) = |(T^a^(r)/(r)e“^‘^(e^ A e^). Then by equating 

these two alternative expressions for the held strength 2-form one gets a set of two 

equations 

dr[a{T)Al{T)] = 0, (10) 

yl“(r)e“'^ + aa{T)A^,{T)Al{T)e^^^ = aa(r)/(r)e“'^^ 


Immediately, one can realize that these two equation can be simultaneously satished provided 

k 


Atir) = 

/M = 


(Ta(r) 
k{k + l) 


( 11 ) 


a^a^r) 

Our next job would then be to demonstrate that this choice of the SO(4)-symmetric ansatz 
for YM gauge potential and its held strength tensor are legitimate since they indeed satisfy 
the Bianchi identity. As a matter of fact, it is straightforward to check that the Bianchi 
identity DF^ = dF'^ — A = 0 is satished with our choice of ansatz, A“ = and 

= ^k{k + l)e“^‘^((T^ A a^). Therefore, we can conclude that a consistent choice of the 
SO(4)-symmetric ansatz for non-vanishing components of YM gauge potential and its held 
strength are 

k 


AKt) = etAl = -Y-X, 

^ aa[T) 


( 12 ) 


_ Hk + l)^abc 


Kir) = ete^F;, = 


a2a2(r) 


Notice, again, that in contrast to the formulation employed in the original work by Hosoya 
and Ogura [5] where they obtained these types of expressions for YM gauge potential and 



its field strength tensor only after explicitly solving the classical YM held equation, here we 
obtained them simply from dehnitions and the Bianchi identity. Therefore, we particularly 
would like to stress that these SO(4)-symmetric ansatz for and remain valid even 
off-shell as well as on-shell. In other words, these expressions for the YM gauge potential 
and its held strength can be used for both classical and quantum treatment that we shall 
discuss later on. As announced in the introduction, the aim of this work is to study the 
physics of classical and quantum wormholes in EYM theory. And we are supposed to study 
the classical wormhole solutions and their ehects on low-energy physics in the present section 
and the quantum wormholes later on. To repeat, wormholes are Euclidean metrics which 
are solutions to the Euclidean classical held equations representing spacetimes consisting 
of two asymptotic Euclidean regions joined by a narrow tube or throat. Thus to study 
classical wormholes, we need to hrst look for exact solutions to the Euclidean classical held 
equations which exhibit above-mentioned geometrical structure. Note that we already have 
chosen SO(4)-symmetric ansatz for the YM gauge potential and its held strength which are 
consistent with the dehnitions and the Bianchi identity. Thus before we attempt to solve 
the Einstein equation, we need to check if these ansatz for YM helds can be solution to 
the classical YM held equation. As we shall see shortly, demanding that they be a solution 
amounts to hxing the parameter “/c” appearing in the ansatz. Now, since we choose to work 
in non-coordinate basis, we also consider the classical YM held equation in non-coordinate 
basis as well. And in non-coordinate basis, the YM held equation is obtained by simply 
replacing F°‘ with its Hodge dual in the Bianchi identity equation given earlier, namely 

DF^ = dP^ - A = 0. (13) 

Since the Hodge dual of the YM held strength is obtained to be 

F“ = A;(A; + l)^^(drAa“), (14) 

a(r) 

it is again straightforward to conhrm that our choice of the SO(4)- symmetric YM gauge 
held given earlier does indeed satisfy the classical YM held equation provided k = —1/2, 
namely if 
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( 15 ) 

= ]-k{k + A a") = A a"). 

2 8 


Next, we turn to the Einstein field equation. Since the S0(4)-symmetric ansatz for YM 
fields are given in terms of the scale factor a(r), the coupled Einstein-YM field equations 
essentially reduce to the Einstein field equation of the scale factor a(r) alone. 

Now consider the rr-component of the Einstein field equation 


Rrr - \g„R + = ^T„. (16) 

The YM field energy-momentum tensor in coordinate basis, T^-, is related to its counterpart 
in non-coordinate basis by Trr = = e°e°Too and Too can be easily computed by 

YM field strength in non-coordinate basis, = 0, = (—to be Tqq = 

Thus also uslug. 


H„=|-3(A+3(Y(-)], 

a JM a 


^= 4 [( 4 )(-) 


a- 


4(-+(-n+4] 

a a 


where overdot denotes the derivative with respect to the Euclidean time r and = 2/37rMp, 


one gets the rr- component of the Einstein equation as 


= (17) 

dr 

where Tq = 371“^/8g‘^ and we introduced “dimensionless cosmological constant” as A = 
16A/9Mp. Now, we attempt to solve this Einstein field equation with different gauge choice 
for the lapse N{t) and in the presence/absence of the cosmological constant. 

1. In the absence of the cosmological constant 
(i) With the gauge choice Y(r) = 1 : 

The Einstein equation above reduces to 


( 




( 18 ) 
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of which the solution is given by 


r = 



r 

Ja{0)=ro 


ada 



which yields upon the integration 


a(r) 




2 ] 1/2 


( 19 ) 


where we chose that the wormhole throat, i.e., the minimum value of a(r) occurs for r = 0, 
namely ^ |r=o= 0 and a(r = 0) = tq. This solution represents a Tolman-type wormhole 
found first by Hosoya and Ogura [5]. 

(ii) With the gauge choice iV(r) = a(r) : 

The Einstein equation, this time, reduces to 



( 20 ) 


of which the solution is given by 


r = / dr' = / 

Jo J a{0)=ro 


which yields upon performing the integration 


da 



a(r) = ro cosh(2r). (21) 

2. In the presence of the cosmological constant 

(i) With the gauge choice iV(r) = 1 : 

The Einstein equation reads 

(^)2 = l_Aa2-d (22) 

dr 

As usual, we choose the minimum value of a(r), i.e., the wormhole neck to occur for r = 0. 
Then from 

|r=o= - r^] = 0 

dr 
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which has zeros at = ^[1 dz yl — 4ArQ], the radius of the wormhole neck, which exists 
when A < l/4rQ, is found to be 


a(0) = r_ 




(23) 


Thus, upon carrying out the integration 



one gets 


= [^{1 - \/l -4Arg cos(2yAr)}]^/l (24) 

This solution for the scale factor a(r) takes the minimum value at r = 0, a^m = [^(1 ~ 
— 4Arg)]^/^ and the maximum value at r = 7r/2V^, a^ax = [^(1 + — 4Aro)]^/^. Thus 

this solution represents a typical wormhole configuration. 

(ii) With the gauge choice iV(r) = a(r) : 

In this gauge, the Einstein equation reads 

^ 2 ^ ^2 _ ^^4 _ 2 _ ^ 25 ) 

dr 


Once again, we choose the minimum value of a(r), i.e., the wormhole neck to occur for 
r = 0, i.e., (^) |r=o= 0, then a(0) = r_ = [^(1 — — dArg)]^'^^ which is possible provided 

A < l/4rQ. Then using (—Aa^ + — r^) = — A(a^ — r^)(a^ — r‘^), 



This integral can be expressed in terms of the elliptic integral and it is 


V^r^r = F{^, ^) (26) 

r_ 

where F{x, k) denotes the elliptic integral of the 1st kind (0 < /c < 1) and F{x = 1, /c) = 
K{k) is the “complete elliptic integral” of the 1st kind. 
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Having constructed wormhole instanton solutions with different gauge choices for the lapse 
N{t), we now turn to their contribution to the topology- changing tunnelling amplitude. 
Namely, we evaluate the wormhole instanton action I eym {instanton) (here Ieym denotes 
the Euclidean action of the EYM theory), i.e., the minimum Euclidean action of the worm- 
hole instanton conhguration which makes dominant contribution to the topology-changing 
tunnelling amplitude. And this amounts to substituting the Einstein held equation (its 
trace) satished by the wormhole instanton solution 

327r 


R = 


M2 

P 


A 


(27) 


into the Euclidean EYM theory action given earlier to obtain 

I EYM {instanton) = f d'^Xy/gl-A + ^{F^^^)'^]. 

J 

Then using, / d'^x^ = {27r‘^a'^) drNa^, = 2/37rMp, A = 16A/9Mp and (E“^)^ = 

= 3/8a'^a‘^, this becomes 


roo 

I EYM {instanton) = / drNa 
Jo 


8A 37r2 1 

+ 


0 

1 /•“ 

2 Jo 


9M4 16c/2a4' 


(28) 


, .N , 4 , 2i 1 fAoo) (^-Xa^ + rl) 
dT{ —)[—Aa -I-rg] = - da ^ 

d Zi J 


—Aa^ -|- 


where we used the rr-component of the Einstein equation {da/dr)'^ = N‘^{1 — Xa? — Tg/a^) 
which yields dT{N/a) = da/ \J —Xa'^ -|- a2 — rg. Note, hrst, that above expression for 
I EYM {instanton) is completely independent of the lapse function N{t). Since the rele¬ 
vant gauge freedom for the case at hand is the arbitrariness in choosing the lapse N{t), 
this means that the wormhole instanton action Ieym {instanton) has gauge-invariance. As 
a matter of fact, since the quantity exp [—Ieym {instanton)] represents the semi-classical 
approximation to the topology-changing tunnelling amplitude, it is a physical observable 
which should have manifest gauge- invariance. And we have just conhrmed this point. 

Next note that unfortunately, it is not possible to carry out explicitly the integral in the 
expression for Ieym {instanton) above to obtain the precise value of the wormhole instan¬ 
ton action. As was pointed out by Hosoya and Ogura [5], however, for a sufficiently small 
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cosmological constant, this instanton action may be expanded in powers of the cosmological 
constant A = 16A/9Mp to yield 

1 Stt 

I EYM (instanton) = —— — -— InfroV^) + 0(V\) (29) 

6 A Aag 

where Tq = 37r‘^/8g‘^ as introduced earlier and ag = The hrst term results from 

the contribution from the half of the de Sitter instanton while the second term represents 
the contribution from the wormhole (neck) geometry especially with the logarithmic factor 
indicating the infrared cut-off generated by the hnite wormhole neck size. We now end 
with few remarks. Firstly, note that the expression for I eym (instanton) given above in 
eq.(29) is valid only for sufficiently small values of the cosmological constant. Next, con¬ 
sider the contribution from the wormhole conhguration to the topology- changing tunnelling 
amplitude 

exp In(royA)]. 

As Hosoya and Ogura [5] remarked, at hrst sight, this factor appears to serve as a suppression 
rather than an enhancement for the tunnelling amplitude since A is assumed to be small 
anyway. Besides, wormholes with larger size, namely, larger than the scale set by the 
dimensionless cosmological constant, i.e., tq >> l/y/X appears to enhance the tunnelling 
amplitude than the ones with smaller size. This last point is particularly disappointing 
since it is just against our expectation. As a matter of fact, it is noteworthy that just the 
opposite happens in the physics of axionic wormholes studied by Giddings and Strominger 
[4] where the tunnelling amplitude is highly suppressed for wormholes with size greater than 
the Planck length. All these arguments might be rather hasty since it is the whole quantity 

exp [—I eym (instanton)] =e^(roV^)^ (30) 

that one really should consider as the right tunnelling amplitude. And for sufficiently small 
A, the hrst singular factor dominates over the second factor (coming from the wormhole 
contribution) still producing a sizable tunnelling amplitude. Neverthless, since it is the role 
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played by the small cosmological constant that essentially elevates the tunnelling amplitude 
to a sizable magnitude, it seems fair to say that the wormhole instanton solutions in this 
EYM theory do not seem to be local minima of the Euclidean action which have manifest 
contribution to the semiclassical tunnelling amplitude. 

III. Quantum Wormholes in Einstein-Yang-Mills Theory 
Thus far, we have considered classical treatment of wormholes in EYM theory. Namely 
we constructed SO(4)-symmetric metric solutions to the classical EYM field equations that 
exhibit the geometry of wormhole spacetimes and studied their nature. We now turn to the 
quantum treatment of wormholes in the same EYM theory. The formulation we shall employ 
can be summarized as follows : we construct and study a minisuperspace model (again 
based on SO(4)-symmetry) of canonical quantum cosmology in which the main objective 
is to solve the Wheeler-DeWitt (WD) equation to hnd the universe wave function. And 
to see if there is an excitation that can be interpreted as a “quantum wormhole”, we look 
for a particular solution to the WD equation with “wormhole boundary condition”. Here 
the wormhole boundary condition refers to an appropriate boundary condition of a universe 
wave function that allows one to naturally identify the universe wave function as representing 
an excitation corresponding to a wormhole state. And as a proposal for such wormhole 
boundary condition, we shall employ the one advocated by Hawking and Page. According 
to their proposal, wormhole wave functions are supposed to behave in such a way that 
they are damped, say, exponentially, for large 3-geometries {y/h oo) and are regular in 

some suitable way when the 3-geometry collapses to zero {y/h 0). We now begin by 

constructing a minisuperspace quantum cosmology model based on EYM theory described 
by the action given earlier. As already stated, to do so, we choose to take the avenue of 
canonical quantum cosmology based on Arnowitt-Deser-Misner (ADM)’s (3-1-1)- space-plus- 
time split formalism [11-13]. As usual, then, in order to render the system tractable, we 
reduce the inhnite-dimensional “superspace” down to a hnite-dimensional minisuperspace 
by assuming that the 4-dim. spacetime has the geometry of spatially-closed {k = -|-1) FRW- 
metric which is the one we adopted in the classical treatment of the system earlier. The 
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geometry of its spatial section is, then, that of and hence it posseses SO(4)-symmetry. 
Now, since the spacetime geometry is taken to possess the SO(4)-symmtry, the matter held 
dehned on it, i.e., the YM gauge held should have the same SO(4)-symmetry. Thus the 
SO(4)-symmetric ansatz for the metric is again the spatially-closed FRW-metric 

ds^ = a^[N\T)dT^ + a2(r)cT“ ® a“] 

= VABe^ ® (31) 


and the SO(4)-symmetric ansatz for the YM gauge potential and its held strength tensor 
are, in non-coordinate basis 


yl“(r) = 0, Atir) = ^-5“ (32) 

aa[T) 

Tpa (^\ — n jpa, ( \ _ fe(fe T 1) 

h, ^ hc\^) 2 2( \ ^ 

where A,B = 0,1,2,3 and a, 6, c = 1,2,3 and again = 2/3nMp. Note hrst that in the 
metric, if N is imaginary, it is the Lorentzian metric with Lorentzian time r = t and if N real, 
it is the Euclidean metric with Euclidean time r. As is well-known and as we shall see later. 


the WD equation and its solution are independent of N and r and hence remain the same in 
both the Lorentzian and Euclidean signatures. Next, recall that the SO(4)-symmetric ansatz 
for A%{t) and given above remain valid oh-shell as well as on-shell as we stressed 

in the classical treatment of the system discussed earlier. And it is because we obtained 


them simply from dehnitions and the Bianchi identity which are bottomline conditions that 


should be met in quantum formulations as well. Thus in the present quantum treatment, 
the parameter k appearing in the ansatz for YM gauge field are undetermined. (As we have 
seen, it gets determined to he k = —1/2 only when one imposes on-shell condition, namely 
only if one demands the classical YM held equation to be satished.) Now, for the sake of 
dehniteness, we choose the Euclidean signature (i.e., Y(r) is real) and write the Euclidean 
action of the EYM theory given earlier in terms of these SO(4)-symmetric ansatz for the 
metric and YM gauge held. Namely, using / d'^Xy/g = JdTN{Jg 3 d^x\/h) = {2T[‘^a^)jdTNa^ 
and 
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( 33 ) 


R 



2 



N 

m 





k\k + iy 

(T‘^a‘^(r) 


with the overdot denoting the derivative with respect to r, the Euclidean action of EYM 
theory takes the form 

- §.R] = i JirNaHA - {i + ^ 

IvM = / = (2-rV‘‘) / dr]Va=>[i{6tY_^}] 


and hence 


— Ig + Iym 




For the sake of completeness we also provide the Lorentzian action of the EYM theory which 
can be obtained via iSeym = —Ieym and r = it 


Seym — -JdtNa^[—\ + {— — = J dtLADM 


where now the overdot denotes the derivative with respect to the Lorentzian time t and again 
A denotes the dimensionless cosmological constant A = 16A/9Mp and Tq = k‘^{k + l)^67r^/5f^ 
which becomes its classical counterpart Tq = 3n‘^/8g‘^ for k = —1/2. Now the subsequent 
procedure toward the derivation of the WD equation via the Legendre transformation to 
the classical Hamiltonian and the Dirac quantization can be best demonstrated in familiar 
Lorentzian signature. Thus from this point on we work with the Lorentzian action for the 
construction of the WD equation. Before we proceed, however, here we make one important 
remark on the structure of the EYM theory action written in SO(4)-symmetric ansatz given 
above. The fact that we choose to consider only the magnetically-charged case = 0) and 
the fact that the SO(4)-symmetric ansatz for YM field strength is written as 


17 



effectively reduce the EYM theory action to the action of a single minisuperspace variable, 
i.e., the scale factor a(r) alone. Thus the resulting WD equation will also become a one-dim. 
Schrodinger-type equation for the universe wave function 'h[a]. As we shall see shortly, this 
reduction to one-dim. problem of minisuperspace model essentially makes it possible to hnd 
the universe wave function which represents a quantum wormhole state. Formally speaking, 
this reduction is achieved by the full gauge-hxing which consists of the gauge choice = 0 
associated with the 3-dim. diffeomorphism in the gravity sector and the gauge choice Aq = 0 
in the YM sector. 

Now getting back to the derivation procedure for the WD equation, we identify the canonical 
momentum conjugate to the minisuperspace variable a with 


QLadm a / .N 

and then obtain the classical Hamiltonian via the Legendre transformation 


Seym — j dtLADM 


(37) 


(38) 


,a — H 


ADM — 


4a - {NHo + NiW)dt] 


where Ni = 0 and Hadm = NHq + NiH^ denotes the ADM Hamiltonian consisting of the 
sum of secondary constraints Hq = 0 and iL* = 0 as we shall discuss shortly. Namely the 
classical Hamiltonian is obtained as 


SSeym 1 3r \ , 1 , 
= xa -A + ^ + 


cS 


5N 


N^a? 


(39) 


= ^[P^-Aa^ + a2-ro^]. 


Namely, 


H. = 


SS, 




where U (a) = — 4 + 


(40) 


General relativity is one of the most well-known constrained system. The invariance of 
the system under the 4-dim. diffeomorphisms (consisting of the time-reparametrization 
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and the 3-dim. general coordinate transformations of the spacelike hypersurface) leads 
to the emergence of 4-constraint equations. Of them, we need not explicitly impose the 
3-momentum constraint equations since we already have taken the iV* = 0 gauge which 
amounts to assuming the 5'0(4)-symmetric spatially-closed FRW-metric. Thus, we only 
need to impose the Hamitonian constraint Hq = 0. The classical Hamiltonian constraint 
now reads 

Hq = + rg)] = 0. (41) 

Next, in order to quantize this EYM system, we need to turn to the “Dirac quantization 
procedure” for the constrained system. According to the Dirac quantization procedure, 
the invariance in the action of the theory under the 4-dim. diffeomorphisms is secured 
by demanding that the physical (universe) wave function T be annihilated by “operator 
versions” of the 4-constraints. Therefore, the classical Hamitonian constraint above turns 
into its quantum version, namely the WD equation given by 

Mva = = 0 - ( 42 ) 

In order to obtain the correct form of this WD equation, we first examine the structure of 
the classical Hamiltonian 

Hq = + rg)] 

= T + V = + V (43) 

where V = ^U{a) = ^(Aa^ — + r^). Here, one can readily read off the “minisuperspace 

metric” Gafs as 

G^p = -a5^f3 , = (44) 

(X 

with 7 “ = a and Hq, = pa being the minisuperspace variable and its conjugate momentum 
respectively. Then, by the usual substitution, 

G^^H^H^ ^ -V" (45) 
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with finally we arrive at the WD equation 


Ho'^ = 


1rl d d , 
2^ada da' 


-17(a) 


'h[al = 0. 


( 46 ) 


Note, here, that the minisuperspace metric (or “DeWitt metric”) Gafsi'y) generally a 
function of minisuperspace variables 7 ". Therefore, in passing from classical to quantum 
version, there arises the “ambiguity in operator ordering” problem. And as have been 
suggested by Hartle and Hawking [12], one may partly encompass this “operator-ordering” 
problem by rewritting 


da da' 


1 d 


d 


aP da^ da' 


(47) 


Finally, the WD equation takes the general form 


-\— — {aF— 
2^aP da^ da' 


+ 17(a) 


T[al = 0 


(48) 


where “p” denotes a suffix representing the ambiguity in operator-ordering. Before we go 
on and attempt to solve this WD equation explicitly, we would like to have some insight 
into the behavior of the solution of this WD equation. To do so, we begin by assigning the 
“normal” sign to the “kinetic” energy term to get 




T[al = 0. 


(49) 


Then the “potential” energy can be identified with 


U (a) = —U (a) = (a^ — Aa"^ — rl). 


(50) 


The Fig. 1(2) given displays the plot of potential energy as it appears in the WD equation 
in the presence (absence) of the cosmological constant A = 16A/9Mp. Note here that 
17(a) = (a^ — Aa"^ — rg) = 0 possesses two positive roots R± provided A < l/4rg or A < 
(|^)^ ^^ 2 (fc+i )2 Afp. Both figures show that due to the contribution to the potential energy, 
{—rl), coming from the YM gauge field sector of the theory, the potential develops a “well” 
in the small-a region for both cases with or without the cosmological constant term. Since 
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the WD equation implies that the total energy of the gravity-matter system is zero, i? = 0, 
the emergence of the potential well in the small-a region reveals the fact that the universe 
wave function \['[a] should be an oscillating function of a there. And this small-a behavior of 
the universe wave function appears to signal the existence of “quantum wormhole” as well as 
other types of spacetime fluctuations in the small-a region of the minisuperspace and hence 
seems consistent with the existence of classical wormhole solution in this EYM theory as 
we have seen in the earlier section. Obviously, the most straightforward way of confirming 
the possible existence of “quantum wormholes” is to solve the WD equation given above for 
the universe wave function which particularly represents an excitation corresponding to a 
wormhole state. Unfortunately, exact, analytic solutions to the WD equation in the presence 
of the cosmological constant are not available. (As a matter of fact, exact solutions to the 
WD equation even for de Sitter spacetime pure gravity are not available either [12,13].) 
However, since our major concern is to hnd the wormhole wave function that is expected to 
exist within the potential well in the small-a region. Erst we focus on the case in which the 
cosmological constant is absent and thus the potential energy takes the behavior depicted 
in Fig.2. Namely we consider the WD equation in the absence of the cosmological constant 



Even in this simpler case, exact, analytic solutions to this WD equation which is valid for the 
whole range of a, [0, cx)) are not available. However, this is not so disappointing since in order 
to see if this WD equation really admits a solution that can be interpreted as a wormhole 
wave function satisfying, say, the Hawking and Page-type wormhole boundary conditions, it 
suffices for us to know just the asymptotic behaviors (such as those at small-a or at large-a) 
of the solution. Thus in the following we consider the behavior of the solution in the small-a 
region and in the large-a region (relative to the scale tq set by the dimensionless coupling 
parameter of the theory, i.e., = k'^{k + /gl). 

(i) For a « tq : 

In this small-a region, the WD equation in the absence of the cosmological constant above 
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can be approximated by 



which admits an analytic solution given by 


'l/[a << ro] 




( 52 ) 


(53) 


Note that this solution in the small-a region is regular at a = 0 provided p < 1 and it 
represents an oscillating solution. 

(ii) For a » tq : 

In this large-a region, the WD equation above reduces to 



which admits an analytic solution given by 


'F[a >> ro] = a*- 2 ^) (55) 

Note that as a —cx), \l'[a >> tq] —thus this solution does involve a damping 
solution in the large-a region. 

Namly the solution to the WD equation in the absence of the cosmological constant pos¬ 
sesses regular, oscillating behavior in the classically-allowed region, a << tq, while it damps 
rapidly enough in the classically-forbidden region, a >> tq in accordance with our expecta¬ 
tion. Then next, in order to see if these asymptotic behaviors of the solution are consistent 
with or equivalent to the “wormhole boundary condition” and hence allows us to identify 
the solution with the wave function of a quantum wormhole, we now refer to a wormhole 
boundary condition advocated by Hawking and Page [10]. Hawking and Page argued in 
a convincing manner that in order for a solution to a WD equation to represent quantum 
wormholes, it should obey certain boundary conditions. And the appropriate boundary con¬ 
ditions for wormhole wave functions seem to be that they are damped, say, expnentially for 
large 3-geometries {-s/h —> cx)) and are regular in some suitable way when the 3-geometry 
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collapses to zero {y/h 0). Particularly, in the context of FRW minisuperspace model, 
large 3-geometries correspond to a —^large limit and the 3-geometry collapsing to zero cor¬ 
responds to a —0 limit. And the damping behavior of the universe wave function at large-a 
indicates that there are no gravitational excitations asymptotically and hence it represents 
asymptotically-Euclidean spacetime while its regularity at a = 0 indicates that it is non¬ 
singular. Therefore, the solution to a WD equation obeying these boundary conditions must 
correspond to wormholes that connect two asymptotically-Euclidean regions. As we already 
have observed, since our solution to the WD equation in the absence of the cosmological 
constant is indeed regular at a = 0 (provided p < 1) and damps sufficiently rapidly at 
large-a, it does satisfy the “wornhole boundary condition”. Namely the WD equation in the 
absence of the cosmological constant does admit wormhole wave function as its solution and 
most probably this will hold true for the WD equation in the presence of the cosmological 
constant. And this is because the essential point that allows the existence of the wormhole 
wave function in the small-a region is the contribution to the potential energy, (—Tq), coming 
from the YM gauge held sector of the theory as we mentioned earlier. Therefore, now we 
can conclude that there are quantum wormholes as well as classical wormhole solutions in 
EYM theory with or without the cosmological constant. 

Next we go back to the case where the cosmological constant is present. We now would 
like to point out that the EYM theory with the cosmological constant may serve as an in¬ 
teresting simple model which provides an overall picture of universe’s entire history from 
the deep quantum domain all the way to the typical, classical domain. To this end, hrst 
note that our knowledge on the nature of the universe wave function in the absence of the 
cosmological constant developed thus far may, in turn, enables us to construct the solution 
to the WD equation in the presence of the cosmological constant at least approximately yet 
quite systematically. Let us go back and consider the WD equation in EYM theory in the 
presence of the cosmological constant given earlier 


1 

2 


\- — — 
- aP da 


“'I;’ + 


Xa^ 



0 . 


(56) 


23 



With the potential energy, U (a) = (a^ — rg), as it appears in the WD equation in the absence 
of the cosmological constant, we now know that the solution to the WD equation represents 
a quantum wormhole states with oscillating behavior in the small-a region. Therefore this 
observation plus the shape of the full potential energy, f/(a) = {a? — \a^ — r^) in the pres¬ 
ence of the cosmological constant as was depicted in Fig.l suggest that the solution to the 
WD equation above would describe the state of the universe that undergoes “spacetime 
fluctuations for very small-a”—“spontaneous nucleation (quantum tunnelling) of the uni¬ 
verse from nothing in a de Sitter geometry” —“subsequent, mainly classical evolution of 
the universe for large-a”. Thus if we are willing to accept {d/da) as the timelike killing 
held in the (mini)superspace, the EYM theory with the non-vanishing cosmological con¬ 
stant appears to serve as a simple yet interesting model which provides a comprehensive 
overall picture of entire universe’s history from the deep quantum domain all the way to 
the essentially classical domain. Then coming back to a practical problem, now we wish 
to construct the approximate solutions to the WD equation. Roughly, the behavior of the 
solution for very small-a will be determined by the wormhole wave function obtained in the 
present work while the behavior for intermediate-to-large-a regions will be approximately 
governed by the de Sitter space universe wave function. Here, the de Sitter space universe 
wave functions have been constructed and extensively studied in the literature associated 
with the issue of initial or boundary conditions for the universe wave function. Therefore in 
order to employ the de Sitter space universe wave function, here we give a brief account of 
the issue of initial or boundary conditions for the universe wave function. The WD equa¬ 
tion is a second order hyperbolic functional differential equation describing the evolution 
of the universe wave function in the superspace. Thus the WD equation, in general, has a 
large number of solutions and in order to have any predictive power, one needs initial or 
boundary conditions to pick out just one solution by, for instance, giving the value of the 
universe wave function at the boundary of the superspace. Thus far, a number of different 
proposals for the law of initial or boundary conditions have been put forward. And among 
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them, “no-boundary proposal” of Hartle and Hawking (HH) [12] and “tunnelling boundary 
condition” due to Vilenkin [13] are the ones which are the most comprehensive and the 
most extensively studied. If stated briefly, the no boundary proposal by HH [12] is based on 
the philosophy that the quantum state of the universe is the closed cosmology’s version of 
“ground state” or “state of minimum excitation” and the wave function of this ground state 
is given by an Euclidean sum-over-histories. Next, Vilenkin’s tunnelling boundary condi¬ 
tion [13] can be best stated in “outgoing modes” formulation which governs the behavior 
of the solutions to the WD equation at boundaries of the superspace. Namely, according 
to this proposal, at “singular” boundaries (such as the region of zero 3-metric and inhnite 
3-curvature {s/h 0)) of superspace, the universe wave function should consist solely of 
outgoing modes carrying flux out of superspace. These two proposals for the law of initial or 
boundary conditions on the universe wave function determine the behaviors of the universe 
wave function like how the universe nucleated (from “nothing”) and then following which 
line it has subsequently evolved to the present one. Indeed in the minisuperspace models 
and within the context of the semiclassical approximation, these two proposals have been 
successfully applied to and tested for simple systems such as de Sitter spacetime in pure 
gravity which is of our interest here or a scalar field theory coupled to gravity. 

Thus employing the well-known de Sitter space universe wave functions with the choice of 
HH’s no boundary proposal or Vilenkin’s tunnelling boundary condition, we now construct 
the approximate solutions to the WD equation. And to do so we consider the WD equation 
in two regions of interest in the minisuperspace. 

Firstly for very small-a region, the WD equation reduces to eq.(51). As we have investigated 
earlier, the important asymptotic behaviors of the solution to this equation are given by 

T/la] —(for a —0) (57) 

—^ (for a large) 

where p < 1 for regularity of the solution. Secondly for intermediate-to-large-a regions, the 
WD equation can be roughly approximated by 
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As mentioned, the semiclassical approximation to the solutions of this de Sitter space WD 
equation has been thoroughly studied. First HH’s no-boundary wave function [12] is given 
by 

exp [i-{l - (1 - Aa^)^^}! (for Aa^ < 1) 

*-3A 

^ (Aa^ << 1), (58) 

_(1±£) r 1 TTn X 2 IN 

= a 2 .lexp [—]2cos -—--— (for Aa > 1) 

OA *- OA 4-1 

_ ( if j ^3 _ j n^~\ /\2 -i\ 

^ a ^ 2 ) e 3 +e ^ (Aa >> 1). 

This HH’s no-boundary wave function consists of both “ingoing” (contracting) and “outgo¬ 
ing” (reexpanding) modes in the classical-allowed region (Aa^ > 1) which, then, decreases 
exponentially as it moves towards smaller values of a in the classically-forbidden region 
(Aa^ < 1). Next, Vilenkin’s tunnelling wave function [13] is given by 

W = ~-r{l ~ (1 ~Aa^ < 1) 

^ OA 

^ (Aa^ << 1), (59) 

= — l)“^/"^exp —4-|l + i(Aa^ — l)^-^^}! (for Aa^ > 1) 

3 A 

This Vilenkin’s tunnelling wave function exponentially decreases as it moves from small 
toward larger values of a (i.e., emerges out of the potential barrier via “quantum tunnelling”) 
and then upon escaping the barrier, it consists solely of “outgoing” (expanding) mode in the 
classically-allowed region (Aa^ > 1). Now we are ready to write down approximate solutions 
to the WD equation in the presence of the cosmological constant by putting these pieces 
altogether. To do so let us denote the smaller and larger roots of the equation 

U (a) = — \a^ — Tq = 0 

by r_ and r+ respectively which exist provided A < l/drg or equivalently A < |[ g^fc(fc_|_i) ]^Afp. 
Then the two roots are given by ^ [1 ± — dAr^]. 

(1) With the choice of HH’s no-boundary wave function : 
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']/[a] = ^/[a] (region / : 0 < a < r+), 

= (region II : < a < oo). 

(2) With the choice of Vilenkin’s tunnelling wave function : 

'l/[a] = \h/[a] (region J : 0 < a < r+), 

= (region JJ : r_ < a < cx)). 

The two kinds of universe wave functions corresponding to the two different choices of the 
boundary conditions are plotted in Fig 3. and 4 respectively. 

IV. Discussions 

Now we summarize the motivation and the results of the present work. 

We revisited, in this work, the Einstein-Yang-Mills theory considered hrst by Hosoya and 
Ogura [5] which is one of the classic systems known to admit classical. Euclidean wormhole 
instanton solution. The classical wormhole instanton in this theory as a solution to the 
classical held equations and much of its effects on low energy physics have been studied 
extensively in the literature. Since this EYM system admits classical wormhole solutions, 
one may wonder if there is any systematic way of exploring the existence and the physics of 
“quantum” wormholes in the same theory. The present work attempted to deal with this 
kind of yet unquestioned issue. In order to explore the quantum wormholes in this sys¬ 
tem systematically, we worked in the context of canonical quantum cosmology and followed 
Hawking and Page [10] to dehne the quantum wormhole as a state or an excitation repre¬ 
sented by a solution to the Wheeler-DeWitt equation satisfying a certain wormhole boundary 
condition. Particularly, in the minisuperspace quantum cosmology model possessing SO(4)- 
symmetry, it is demonstrated that there indeed exists a solution to the Wheeler-DeWitt 
equation whose asymptotic behaviors satisfy the appropriate wormhole boundary condition 
for the case when the cosmological constant is absent. Thus we confirmed our expectation 
that the EYM system does admit quantum wormholes as well as classical wormholes. Fur¬ 
ther, we pointed out that the minisuperspace quantum cosmology model based on this EYM 
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theory in the presence of the cosmological constant may serve as an simple yet interesting 
system displaying an overall picture of entire universe’s history from the deep quantum do¬ 
main all the way to the classical domain. 

As we have stressed in the text, the essential point that allowed us to explore, in a concrete 
manner, the quantum wormhole in the context of the minisuperspace quantum cosmology 
model was the following observation. In their original work, Hosoya and Ogura [5] ob¬ 
tained the expressions for the YM gauge potential and its held strength tensor only after 
explicitly solving the classical YM held equation. In particular, the SO(4)-symmetric YM 
held strength tensor obtained in this way, = 0, in turn, reduces the 

Einstein equation to that of the scale factor a(r) alone. However, we realized in this work 
that even without imposing the on-shell condition (i.e., the classical held equation), one can 
“derive” dehnitions and the Bianchi identity. Therefore this 

SO(4)-symmetric ansatz for the YM held strength remains valid even oh-shell as well as 
on-shell and hence can be used in the quantum treatment of the EYM system. Consequently, 
the Wheeler-DeWitt equation in the context of the canonical quantum cosmology becomes 
a Schrodinger-type equation of the minisuperspace variable a (the scale factor) alone and 
can in principle be solved particularly in the absence of the cosmological constant. 

Finally, our study of quantum wormholes in this work appears to demonstrate that, after 
all, the EYM theory is a simple yet fruitful system which serves as an arena in which we 
can envisage quite a few exciting aspects of quantum gravitational phenomena. 
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Fig.l The plot of the potential energy as it appears in the WD equation in the 
presence of the cosmological constant. 



Fig.2 The plot of the potential energy as it appears in the WD equation in the 
absence of the cosmological constant. 








Fig.3 A regular solution to the WD equation in the presence of the cosmological 
constant employing HH's no-boundary proposal. 



Fig.4 A regular solution to the WD equation in the presence of the cosmological 
constant employing Vilenkin's tunnelling boundary condition. 









